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1. Introduction 

Quantum particles on the real line interacting through a 5-potential are governed by 
the Hamiltonian 

j=i i j<k 

The number of particles, n, is fixed throughout and x = {xi, . . . ,Xn) denotes the 
positions of the particles. We will restrict ourselves to the bosonic subspace of symmetric 
wave functions. Eq. f ll.ip is the Hamiltonian of a quantum many-body system which 
can be analyzed through the Bethe ansatz. The repulsive interaction, k < 0, has been 
studied in great detail and we refer to |ll[2l[3l|4l|5l|6l|7]. The attractive case, k > 0, has 
received less attention. One reason is that the structure of the Bethe equations is more 
complicated. On top, physical applications are not obviously in reach. In the recent 
years, there has been renewed interest. We have now available a detailed study of the 
eigenfunctions [H El [TOl |TT] and, as argued by Calabrese and Caux [12], applications 
to real materials are in sight. A further motivation comes from the one-dimensional 
Kadar-Parisi-Zhang (KPZ) equation [13]. Its replica solution is given in terms of the 
propagator of the attractive 5-Bose gas [lU [15] which can be used to obtain exact 
solutions for some special initial conditions [l6l[I71[l8l[IIl[T9l|2ni[ail22l[23l[2ai25]. 

In the KPZ context, and also in other cases, one is actually interested in the 
quantum propagator (x|e~*^'' |?/), t > 0. In principle, e"*^*" can be expanded in a sum 
(integral) over eigenfunctions. But one might hope to have at disposal more concise 
expressions for the propagator. In the repulsive case, Tracy and Widom [26] carried 
out such a program. The resulting expression we refer to as TW formula, which will 
be discussed below, including its relation to the expansion in eigenfunctions. A natural 
issue is to extend such a program to the attractive case, which is the topic of our 
contribution. 

By symmetry the propagator (x|e~*'^'^ |?/) can be restricted to the domain A = 
{x\xi < . . . < Xn} C M". Using the Bose symmetry, if^ of (11. ip is then defined by 

H^i:{xi,...,Xn) = - '^-^ip{xi,...,Xn) , X E A° , (1.2) 



with the boundary conditions 
d d 



dxj+i dxj 



+ k] ij{xi,...,Xn)\^,^^^^, = , (1.3) 



where the limit taken from the interior, A°, of A. The Hamiltonian is 

a self-adjoint operator and (xle'^^^ ly) is continuous in x,y E A. In particular, 

1 " 

lim(a;|e-*^'= 1^) = ^ IT ^^^i " ^i) ' 2/ ^ A . (1.4) 

i=i 

Throughout the paper x,y G A, hence the position of the particles are ordered 
increasingly. As will be proved in Appendix A, (x|e~*^'= ||/) is analytic in k for otherwise 
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fixed arguments. Thus our strategy will be to analytically extend the TW formula, valid 
for K < 0, to K > 0. As written, the TW formula becomes singular at k = 0. Therefore 
the main task is to understand the structure of the analytic continuation in n. As a 
result, we will arrive at various formulas for the propagator. One formula will be just 
the expansion in Bethe ansatz eigenfunctions, which thus implies their completeness. 

The issue of completeness for the attractive (5-Bose gas on the line has been 
studied before. In his thesis, Stephen Oxford [27j proves completeness of the 
generalized eigenfunctions defined as bounded Bethe ansatz eigenfunctions. He uses 
functional analytic methods to construct the Hilbert space isometry from the generalized 
eigenfunctions and thereby the spectral representation of if^. A similar strategy is used 
by Babbitt and Thomas [28] for the ground state representation of the ferromagnetic 
Heisenberg model on Z. Heckman and Opdam p9] exploit the fact that the 5-Bose gas 
turns up in the representation theory of graded Hecke algebras. (We are grateful to 
Balazs Pozsgay for pointing out this reference.) They have results for the case when the 
interaction strength is allowed to be pair dependent. But only for their expression 
simplifies and they arrive at a Plancherel formula, which is the completeness relation. 

For the system on the line, studied here, the set of admissible wave numbers 
is known explicitly. For a bounded system, in particular with periodic boundary 
conditions, the discrete set of wave numbers are the solutions to the Bethe equations, 
a coupled system of n transcendental equations. Completeness becomes more difficult 
to establish and to our knowledge only for the repulsive case a completeness proof is 
available [30l| . 

The article is organized as follows. In Section [2], we recall the Tracy and Widom 
formula for the propagator in the repulsive case, and rewrite it in terms of Bethe 
eigenstates. In Section [3l we summarize our main results on the propagator with 
attractive interactions. These results are proved in Section HI by performing explicitly 
the analytic continuation to k > 0. A further rewriting represents the propagator in 
terms of the known Bethe eigenstates. The special case of the propagator with all 



particles starting and ending at is handled in Section [5l In Appendix A, we prove 



that the (imaginary time) propagator is an analytic function of the coupling. 
2. (5-Bose gas with repulsive interaction {k < 0) 

Let Sn be the set of all n! permutations of the integers between 1 and n. In the following, 
we use the notations 

n- n (2 1) 

j<k l<i<fc<n 

For K < the TW formula states (in the notation of [26] the strength of the potential 
is called c = — /t) 

<^l^"''''l^^ = d(^I/^---^^" 



The propagator of the attractive delta-Bose gas in one dimension 4 

n . n 



<T(j) > <T(fc) 

The connection to the eigenfunction expansion can be seen by symmetrizing over all 
permutations of the qj. We introduce a new permutation r G and replace qj by qrij)- 
Replacing then a by r^^ o a, one finds 

(x|e-*^'=|?/) = ^ V / dgi . . . dg. 



7t . /t 

n + n (e^^^(--i(.)-^-^0))e-*''.') . (2.3) 

, , 5'? ~ ~ . -I ^ ^ 



One notes the factorization 

A - gfc + i/t _ -A- qj - q^ + in -A- - qk - in 

■^"'(j) > <^"'(fc) <T-l(j) > £r-l(fc) £r-l(j) < <T-l(fc) 



?1 



We introduce the Bethe eigenstates of the Hamiltonian ( ll.ip with repulsive interaction 
with momenta q = (gi, . . . , g„) G M". The eigenstate ipi^'i <?) has energy 

n 

E{q) = Y.^l (2.6) 

In terms of Bethe eigenstates, Eq. (12. 3p . (12. 4p rewrite as 

(x|e-*^''|y) = / dgi ... dg„ 7^(2;; g)^<^e-*^(^) , (2.7) 



where (. . .) denotes complex conjugation. At t = 0, (12.71) reduces to the completeness 
relation for the Bethe eigenstates, 

^ = 7^ / ...dg„|V^(g))(^(g)|. (2.8) 

A proof of the orthonormality of the Bethe eigenstates can be found e.g. in [TT] , 
Appendix A. 
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3. 5-Bose gas with attractive interaction (/t > 0) 



As shown in [Appendix A[ the propagator is an analytic function of the coupling k for 
t > 0. By analytic continuation of fl2.2p from k < to k > 0, we will derive in Section 
[3] an exact expression for the propagator in the attractive case k > 0. 

Before stating the two main theorems, a few definitions are needed. We call 
Dn^M the set of the M-tuples n = (ni, . . . , um) such that Uj > 1, j = 1, . . . , M and 



ni + 



+ um = n. For n G DnM, the clusters ^j{n) 



, . . . ,M are defined by 

uj — {ni + . . . + nj_i + 1, . . . ,ni + . . . + Uj} . (3.1) 

From the Bethe ansatz point of view, the clusters will correspond to bound states of the 
particles. The function = r, acting on {1, . . . , n}, is defined by 

r(a) = s for a = Ui + . . . + nj_i + s E Qj . (3.2) 

More visually, one has 



a 


1 ... rii 


rii + 1 


... ni + n2 




n — nM + 1 • • • n 


Cluster 












r(a) 


1 ... rii 


1 






1 ... Um 



(3.3) 



Finally, we call S'^{n) (respectively S'^i^n)) the subset of Sn containing only the 
permutations a (resp. r) such that for all j = 1,...,M and a,b & Qj with a < b 
one has cr~^(a) < (J~^{b) (resp. T~^{a) > T~^{b)). 

In the attractive case, the following expression for the propagator is proved in 
Section HI 

Theorem 1. For fixed n, let fij , j = 1, . . . , M , he arbitrary real numbers satisfying the 
constraint 



— Uj < /ij < . 
For K > and x,y & A, one has 

n n-M ^ r 

M=l ^ neD„,Mi=l 

M 

X] n n (^e'^*"^'''^^^'^''^"^""^^^^^'<'-'(°'"^--'w^e"*^''^+'''^''J+''^"^^ 

ae5;{n) r€Sii{n) j=l aen, 
M 

X n n 



(3.4) 



j, k — 1 aSSlj 
j k b e Qf. 

<T-l(a) > <t"1(6) 



{qj + 'm{(jij + r(a))) - (g^ + 'm{^ik + r{b))) + in 
{qj + ininj + r(a))) - (g^ + ifi:(/ifc + r(6))) - in 



. (3.5) 



yl// i/ie apparent poles in the integrand cancel except for simple poles at qj + inHj = 
qk + iK^fik + n/c) and qj + 'm{jj,j + rij) = qu + infJ'k, j < k. The integrand vanishes at 
qj + iKfij = qk + ii^fJ^k cind qj + in^iij + Uj) = qk + ii^ifJ^k + ^fc)- 
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Compared to the TW formula (I2.2p . our result (13.51) is more complicated. Physically, 
the complication can be traced to the presence of bound states for the attractive case: 
the n particles are arranged in M clusters of size ni , . . . , um , hence the extra summations 
over M and n. Furthermore, Eq. (13. 5p contains a summation over two permutations 
a and r instead of only one for the TW formula (12.21) . In the special case x = y = 
discussed in Section [5l both summations over a and r can be eliminated. 

For the attractive case, the propagator can also be written in terms of a summation 
over the eigenstates of the Hamiltonian (II. ip . The Bethe eigenfunctions for attractive 
interaction are (see [Hj, Eq. (B.26) and (B.48); in [11], r(a) is equal to r{a{a)) with 
our notations, and Qj to a^^{Qj)) 

n-M M M 



fr i-r / [q, + i/c(r(a) - ^)) - {q, + i^{r{b) - f )) + 

<T-l(a) > <T-l(b) 

with M = l,...,n, n e L'„,M and g G M^^. Eq. ([MD is an eigenfunction of the 
Hamiltonian (II. ip with eigenvalue 

M / 2 
2 



E{M,n,q) = Y,l 
i=i ^ 



njQj - Y^iUj -Uj) ) . (3.7) 

The relation of the propagator with the Bethe eigenfunctions is stated as next theorem 
in terms of (13. 6 P and (13. 7p . 

Theorem 2. For k > and x,y & A, one has 

" 1 f 

M=l ^ ^ neZ)„,M"^*^ 

^(x; M, n, g) ^{y;M,n, q) e-*^(*^'"'«) . (3.8) 

As in the case of repulsive interaction discussed in Section |2l taking t = yields the 
completeness relation for the Bethe eigenstates (13. 6p . Their orthonormality is proved in 
dU, Appendix B. 



4. Analytic continuation from k < to k > 

In this section, the TW formula (12. 2 p for the propagator is extended by analytic 
continuation to the attractive case k > 0. Theorem [T] and Theorem [2] are proved. 



Contribution of the residues 

The contours of integration in (12.20 can be moved freely as long as the denominators 
Qj ~ Qk — keep a strictly positive imaginary part. In particular, if the integration is 
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shifted to e IR + iA(n — j), j' = 1, . . . , n, with A > 0, we obtain a formula vahd for all 
K such that '^{qj — Qk — i'^) — {k — j)X — k > for j < /c, i.e. for all k < X. One obtains 

n . n 

^ n Ff^n(«""'^-'"-''"'«-"")- (^-^ 

crGbn a < b ^ ^ a=l 

In the following, we want to further move the contours of integration, but this time 
the contours will have to cross poles of the integrand, which will add several new terms 
resulting from the residues at these poles, symbolically. 




(4.2) 



If one denotes by j — )■ the action of taking the residue at qj — Qk + i^'^ with r e Z, 
the terms, obtained after moving the contours of integration, correspond to collections 
of j — 7> /c such that, for each £=l,...,n, £—)■... appears only once in the collection 
(since after taking the residue at qi = Qm + ^i^f^, the integrand no longer contains qg). 
Each term thus corresponds to a forest (a set of trees), for example 

1 
I 

{1^2,2^7,3^5,4^7} ^ 2 4 3 (4.3) 

\ / I 

7 5 6 

The particular trees obtained in this fashion depend on the order in which the contours 
are moved. Here, we choose to move first the contour for qn-i in such a way that it 
crosses only the pole at = qn + i/«- Then, we move the contour for g„_2 in such a 
way that it crosses only the poles at g„_2 = qn + ^i^ and g„_2 = ^n-i + (in which case 
we still have an integration over both and qn), or only the pole at g„_2 = qn + 2ifi: 
(in which case the residue at qn-i — qn + ^K, has been taken). We continue in this fashion 
until in the final step the contour for gi is moved. 

In principle, after moving the contours, the propagator (x|e~*^'" \y) will be expressed 
as a sum over forests. In fact, it turns out that during this procedure there are many 
cancellations which remove all the forests which contain trees with "branches" : in other 
words, only the forests with merely "branchless" trees (like a, a ^ b, a ^ b ^ c, 
a ^ b c ^ d, . . .) remain after these cancellations. Instead of a sum over 
forests, we end up with a sum over partitions of {1, . . . , n} (each element of the partition 
corresponding to one of the branchless trees of the forest) . 

In the context of the distribution of the leftmost particle in the asymmetric 
simple exclusion process, the procedure described here bears some similarity with the 
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transformation from Theorem 3.1 to Theorem 3.2 in [31j, where contours of integration 
are moved from small to large circles. A complication in our context is that we need a 
summation over all partitions of {1, ... , n} and not just over subsets of {1, ... , n}. We 
expect that in the case of the full transition probability for the asymmetric exclusion 
process an expression with an integration over large circles would require summing over 
all partitions of {1, . . . , ra}. 

A proof of the previous statements is based on induction w.r.t. an integer I such 
that all the contours for g^+i, . . . , g„_i have already been moved. 

We introduce a few notations. For a boolean condition c, l{c} is defined to be 
equal to 1 if c is true and otherwise. For n G Dn^M, the set Pn{n) contains all the 
partitions A = {Ai, . . . , Am} of {1, . . . , n} with \Aj\ = rij, j = 1, . . . , M. The partition 
A verifies U . . . U A]\j = {1, . . . , n} and for j ^ k Aj fl = 0. The partitions are not 
ordered, i.e. the partition B = {Ajk^i^, . . . , Aj^^m)} is considered to be the same element 
of Pn{n) as A for all R G Sm- Each Aj is called a cluster, and will correspond to a 
bound state of particles in the Bethe ansatz point of view. For a partition A, we define 
d^{a), a = 1, . . . , n, (abbreviated as d{a) to lighten the notation) to be the rank of a in 
its cluster Aj, starting with rank for the largest element of the cluster, rank 1 for the 
second largest, . . . , and rank — 1 for the smallest element of Aj. 

With these notation, the following lemma can be stated. 

Lemma 1. Let i be an integer between and n — 1. For fixed M = 1, . . . ,n, let ej, 
j = i + 1, . . . , M be distinct numbers with < < 1. Then, for < k, < X one has 

M E nK'K-i)!)n(/ 




neD„^M i=i i=i 
e M 

Y Il'^{Aj=m n n M'^-'H>^-'m 

A€P„(n) <^^S„ j=l j=l a,b e Aj 

a < b 



M 



j=l aeAj 
M 

X 



nTT / ilj + i^^d{a)) - {qu + iKd{h)) + in \ 
J-i \ {qj+iKd{a))-{qk + iKd{b))-iK) ' 



J, fc — 1 a € Aj 
j T^k be A^ 
a < b 



Proof. The constraint on the e^-, j = £ + 1, M, implies that (14. 4p is well defined since all 
the poles are at qj = qk + i5 with ^{qj) ^ '^{qk) + 

For i = n — 1, the identity between the expressions (14. ip and (14. 4p of {x\e~'^^''\y) 
is immediate. All the clusters must have size 1, and only M = n contributes. Since the 
poles for g„ are at g„ = qj — ik, j = 1, . . . ,n — 1, the contour for g„ can be moved freely 
from M to M — iKe,„, provided < < 1 and k, > 0. 

We now proceed to prove the general identity by induction in i: we assume that 
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the expression ( 14. 4 p for {x\e~*^'^\y) is valid for £ > 1 and will establish the expression 
with ^ replaced by £ — 1. 

For given n, we want to move the contour of integration for qi from M + iA(n — £) to 
M — i/te^ with < < 1 and different from all the other e^, A; = £ + l,..., M. In order 
to accomplish this, one needs to take into account the residues of the poles dX qi = z for 
—nei < ^{z) < X{n — i). The only poles for qi are at z = qj — in, j = 1, . . . , i — 1, and 
aX, z = qm + ii^{d{c) + 1), m = £ + 1, . . . , M, c G A^. In the first case, using k < A and 
i < ^ — 1, one finds '^{z) > \{n — £), which implies that these poles do not contribute 
when moving the contour for q^. In the second case, using 0<k<A, 0<em<l and 
i + d{c) + 1 < n, one has — /te^ < < '<S{z) < X{n — £), which implies that all these 
poles contribute a residue (with a factor — 2i7r corresponding to a clockwise contour 
integration) . 

Moving the contour for qi produces several terms: one term corresponding to the 
integration over G M — i/te^, for which the integrand still depends on qi, and one term 
for each c G Am, m = £ + 1, . . . , M, for which the residue at qt = qm + iK{d{c) + 1) 
has been taken. The latter term corresponds to merging the cluster Ai = {£} and the 
cluster Am- Assuming > cr^^(c) (otherwise, the pole vanishes), this term is equal 

to 

(-2i^);;T7^ n(^^'(^. - 1)0 n / dg, n / ^g, 



M 



^gP„(n) '^£•5'" i=l i=l b 6 

a < 6 



{a-i(a)>a-i(fe)} 



gi(g„i+iK(d(c)+l)){a;^_l|.^j-j/f) ^_t(g„+iK{d(c)+l))^ 



n n 



^i(qj+iKd{a))(x^.X(^-^-ya)^-t{qj+iK.d{a)f 



] = 1 



M 



X 



n n 

J, fc — 1 a £ Aj 

J # fc i) G Afc 

J, fe / « a < b 

<T-l(a) > <T-l{i,) 

A/ 

n n 

j,k — 1 a G Aj 

j ^ k a < . 



j 7^ « <T-l(a) > <T-1(«) 
A/ 



X 



n n 

, fe = 1 6 e 

J 7^ fe £ < 6 



(gj + iKd{a)) - (g^ + iKd{b)) + Ik 
(g^ + i/t(i(a)) — (g^ + iKd{b)) — in 



{qj + iKd{a)) - {qm + iK{d{c) + 1)) + Ik 
(gj + i/t(i(a)) — (g^ + ii^{d{c) + 1)) — i/t 

(gm, + iK{d{c) + 1)) - (gfc + iKd{h)) + i/t 
(gm + \K{d{c) + 1)) - (gfc + 'md{h)) - m 
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{qra + m{d{c) + 1)) - {qm + md{h)) + in 



X 



n 



b g A„ 
6 7^ c 
1(0 > £r-l(b) 



(g™ + iK((i(c) + 1)) - {q„, + - i/t 



X + iK{d{c) + 1)) - (g^ + iKd{c)) + iK) . (4.5) 
The last line of (14. 5 p contributes a factor 2i/t and the line before contributes 

f djc) - d{b) + 2 \ 

11 ( dic)-d{b) ) • 



6 e Am 

b ^ c 

-He) > <y-Hb) 



Let us first assume that c = min(Am). Then, for all b G A^, b ^ c one has 
(T^^(6) < (T^^(c). Together with cr~^(c) < (T^^(£), it implies that all the elements of 
the cluster (except c) contribute in (14. 6p . This results in a factor (n^ + 1)^^/2. 
Combined with (— 2i7r) and 2iK, we obtain a factor 2TTK{nk + l)^fc. The term with 
c = min(y4m) thus corresponds exactly to the term of (14.40 with £ replaced by £ — 1 and 
the partition A replaced by obtained from A by merging the cluster {£} with Am 
(after a renaming of the gj, ^j, to g^-i, Cj-i for £ + 1 < j < M). 

It remains to show that for c 7^ min(74m), the residues cancel each other. Since 
the cr^^ib) are ordered in the same way as the d{b) for b G A^, there exists a unique 
number / G such that for b G A^, if & > / then cr^^{b) < cr^^(£), and if 6 < / 
then cr~^(6) > cr^^(£). Since cr~^(c) < (j^^(£), one has necessarily / < c (or equivalently 
d{f) > d{c)). Then, (14. 6 p rewrites 

' d{c) - dip) + 2^ 



b G A„ 
ij 5^ c 



The rest of the argument depends on the relative values of d{c) and d{f). If d{f) > 
d{c) + 2, then, there exists b G A^ such that b > f and (i(6) = (i(c) + 2, thus (14. 7p is 
equal to zero. Since d{f) > d{c), the only cases left are d(f) = d{c) + 1 and d{f ) = d{c), 
for which (14. 7p rewrites respectively 

d(c)-dif) + 2 ^ f djc) - djb) + 2 \ (^(c) + l)(rf(c) + 2) 



d{c)-d{f) ^11 V rf(c)-t/(6) 



and 



yr f d{c) - d{b) + 2 \ _ {d{c) + l){d{c) + 2) 

11 1 rf(c)-rf(6) )- ^ ■ ^^-^^ 



6 e A„ 

b > c 



Let us call c' the element of Am such that d{c') = d{c) + 1 (c' is the smallest element of 
Am larger that c). One notes that the two previous cases are exchanged when replacing 
a^^ by o 6i c', with 6'^ the permutation exchanging i and c'. Thus, summing over 
all permutations a, the residues at qt = qm + ii^{d{c) + 1) cancel. □ 
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Our construction achieves the proof of fl4.4p for < k < A and i between and 
n — \. In particular, for £ = 0, one has 

" n-M ^ 

(x|e--1,) = 5: ^^-^ E riK'K -l)!) 

M=l ' neD„^Mj=^ 

M / „ \ M 



X n / z z n n ^w-Ha)>a-Hb)} 

a < b 

M 



j=l aeA 
M 

X 



nTT / (qj + indja)) - {qt + iKd{h)) + i/t ^ 
\{qj+i^d{a))-{q, + iKd{h))~iK) ■ ^ 

'j be aI 

a < b 

This expression no longer depends on A. Hence, it is valid in the entire range k > 0. 
4-2. Partitions and permutations 

Exchanging Aj and Ak in f l4.10p is the same as exchanging qj and qk, or ej and e^. Since 
the ej are arbitrary numbers satisfying a constraint (0 < e^- < 1 and all ej different) 
which is the same for all j, it is possible to add an extra sum over all permutations of 
the Aj, compensated by a global factor 1/M!. This is equivalent to summing now over 
ordered partitions A = {Ai, . . . , An), such that for all R G Sm different from the identity 
permutation, the ordered partition B = (v4/j(i), . . . , Ar(^m)) is distinct from A. 

There exists a bijection between ordered partitions A such that \Aj\ = Uj, 
j = 1,...,M, and permutations r G S"{n) (5'^'(^) is defined after Eq. (13. 5p ). By 
this bijection, the cluster Aj is equal to {r~^(a),a G Qj{n)} = T~^{Qj) and one has 
1 + d^{a) = r^(r(a)), using the definitions (13. ip and (13. 2p . Eq. (I4.10p becomes 

" n-M ^ 

M=l ^ ' n€D„^M j='^ 

M / „ \ M 

X n / z z n n ^w-Ha)>a-Hb)} 



a < b 

M 

i(gj+iK(r(r(a))-l))(x^_i(^)-S/a) -i(qj+iK(r{T{a))-l))= 



j=l agr-i{nj) 
M 

X 



TT TT / {qj + i/tr(r(a))) - (g^ + i/tr(r(&))) + ift: \ 

^^},A{Qj + ^^r{r{a)))-{q, + iKr{Tm-i^J ■ ^ 

J ! — -L a t T 'j / 
JT^fc i>Sr-l(nfc) 
a < 6 
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One can now replace everywhere a and h by T^^(a) and r"^(6). We also replace a by 
OCT. Because of the definition of S'^{n)^ if a, 6 G VLj with T^^(a) < t^^(6) then a > b. 
This implies that the constraint with the is equivalent to cr G S'^{n) (defined after 
Eq. (I33D). We obtain 

n-M M f r 

M=l ^ ^ ^6D„_„j=l j=l \^M-iK£, 

M 

5Z 5Z n n (e'^*+'''^''^"^"^^^^''-''w"^--'(°)^e"*(''^+'''(''^")"^^) 

(7G5;{n) TeS;'{n) i=l aGHj 
X 



(xle-*^' 



J, fc — 1 a £ Hj 



^.c?. Po/e structure of the integrand and summation over Bethe eigenstates 

{qj + i/tr(a)) - (g^ + iKr{b)) + in' 



One has the factorization 

M 

n n 



j, k — 1 a £ f2j 

"1(a) > <T 
-1(a) < r-l(6) 



(^j + iKr{a)) - {qk + iKr{b)) - Ik 



n n 

j<k a e fij 

<T-l(a) > <T-l(b) 
T-l(a) < T-l(6) 



(gj + iKr(a)) - (g^ + iKr{b)) + ik 
{qj + i/€r(a)) - (g^ + iKr(6)) - in 



X 



^(gj + i/€r(a)) - (gfe + i/tr(6)) - i/t 



n n 



j<k a e sij 
b e sik 

<T-l(a) < <T-l(b) 
T-l(a) > T-l(b) 



(gj + i/tr(a)) - (g^ + iKr(6)) + 



IK 



Y[ Yl ^ ^ i«r(a)) - (gfc + iKr(6)) + i/t 

<T-l(a) > (T-l(b) 



(gj + iKr{a)) - (g^ + mr{h)) - in 



M 

X 



j<fe a e n, \{(lj + '^>^r{a))-{qk + iKr{b)) + iK 

b e Ofc 



We introduce the functions 

n-M M M 



nTT / (Qj + i^r(a)) - (gfc + i/tr(&)) + ift: \ 
J-J- \ {qj + iKr{a)) - {qk + iKr{b)) - in I ' i - ) 



nTT / (Qj + i^r(a)) - (g^ + iKrjb)) - ift: \ 
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n*^ -TT / {qj + mr{a)) - (g^ + mrjh)) + iK \ 



X 



6 G fife 

and 

n-M M M 

i=l reS;'(n) i=l aGCj 

^ TT TT (gj + ift:r(a)) - {qk + i/tr(fe)) - m \ 

One notes that r G 5'"(n) is equivalent to R{n) o r G with R{n) the permutation 

that exchanges a, 6 G j = 1, . . . , M, if r(a) + r(6) = Uj. This imphes 



M, n, g) = M,n,q- inn + ik) , (4.16) 



(x|e-*^« 



where (...) denotes complex conjugation. Thus, one can write 

M 

II II (^e-'^i.+-Mria)-i)?^ n, q)<f,{y;M,n,q-i^n + i^) . (4.17) 

j=l a£Qj 

At this point, one wants to move again the contours of integration so that the exact 
Bethe eigenfunctions appear. One then needs to know precisely the location of the poles 
of the integrand in fl4.17p . To achieve this, we prove the following lemma: 

Lemma 2. The expression 

^.(x;M.g.,)nnn( ^'-^-+y7'-/> ) . (4,i8) 

]<kr=ls=l ^ -"J ^ ' ' 

with iPk{x; M,n,q) defined in Eg. (^77^, is holomorphic as function of the qj 's in the 
whole complex plane. 

Proof. One defines xi^'i M,n,q), equal to f l4.18p up to a global normalization 
independent of g, by 

M 

x{x;M,n,q)= JJ JJ (^e^fe+^W-i^K-K.)) 

TT TT / (gj + iKr(a)) - (gfc + iKr{b)) + i/tsign(ff"^(a) - a~^(6)) \ , . 
}i V (g,+iKr(a))-(g, + i«:r(6)) J - ^ - ) 

b G Ofc 

Noting that the product 

TT f ilj + i/tr(a)) - {qj + i/tr(&)) + ift:sign(a-^(a) - a~^(&)) \ ,^ ^q) 
V (g,+i/€r(a))-(g,+i«:r(6)) J 

a < b 
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is nonzero for all j = 1, . . . , M if and only if the permutation a is an element of S'^{n), 
one can instead sum over all permutations a & Sn- One has 

M M 

x{x;M,n,q)=l[n,\ (e^^^+^'^'-^^^-^'^^^-^w) 

i=l (TG5n j=l aeClj 

M n 

X 

j<k a<b 

Then, one defines x(,x;M, n, ^) by replacing everywhere in M, n, g) the variables 
Qj + ift;(r(a) — 1), j = 1, . . . , M, a e flj by with the result 

M 



AAA-L^ (g,+i«:r(a))-(gfc + i«:r(6)) /' ^' ^ 



X 



-Q / Ca-6 + i/^sign((7 ^(a)-a ^(&)) \ 

a<b ^ - 6 / 

If x(a;;M, n, ^) is analytic in the ^a, then x(,x;M, n, will also be analytic in the ^j. 
One notes that x(a;; M, n, ^) has only simple poles. The term with the permutation 
a & Sn oi the residue at — with c < d is given by 

M 



(ivi \ n 

j=l ^ a = 1 



X 



X 



n 

a <b 
a ^ c,d 

n 

n 

6=1 



n 

X 



6=1 
b^c,d 











c) - a-\d)) 






- ib + IK sign{a-^ (a) 






ia - ib 




id 


- ^6 + i«;sign(cr~^(c) 






id - ib 




id 


- ib + msign{a~\d) 




id - ib 



This term cancels with the term for which a~^{a) and (J~^(&) are exchanged. Thus, 
x{x; M,n,i) and x(x;M, n, are analytic functions in the ia and in the qj. This 
finishes the proof of the lemma. □ 

Then, since 

-Q -Q /^gj - gfc + i/t(r - s - 1) Qj - qk + iniuj - - r + s + 1) 

r=l s=l 



Qj -qk + ii^[r - s) Qj -qk + iK[nj - rik - r + s) 

(4.24) 



(Qj -qk- iKnfc)(gj -qk + ii^rij) 



{Qj - Qk){qj -qk + ii^irij - Uk)) ' 
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the integrand in (14.171) is analytic in the qj except for simple poles at qj = qk + inuk 
and qk = qj + i^rij for 1 < j < k < M. At this point, it is possible to set ej = 
for all j = 1, . . . , M in f l4.17p since the poles at qj = qk cancel. Then, the contours 
of integration can be moved again, from qj G M to G M + i/ij. Before moving the 
contours, one has (with k, > and j < k) 

'^{qj — qk — ii^nk) = —KUk < and '^{qk — qj — i^Uj) = —urij < . (4.25) 

After moving the contours, these inequalities must still be valid. One notes that if the 
/ij's satisfy the constraint 

- rij < fij <0 , (4.26) 

both inequalities are satisfied. Shifting the contours and making then the change of 
variables qj — )■ qj + infij in (14.121) leads to the result (13.51) of Theorem [1], while shifting 
the contours and making the change of variables in (14.171) leads to 

n „ M 

M=l ^ ^ n&D„,M i=laGn,- 



g-t(g,+iK(Mj+r(a)-l))' 



X (Pk{x; M,n,q + mjl) (p^iy, M, n,q + injl + inn — in) . (4.27) 
One notes that 

^ {qj + + r(a) - 1))^ = + Imq^ ^(/ij + r(a) - 1) - + r - 1)^ 

aGf2j r=l r=l 

= n^q] + 1mq,n, + ^""^'"^^ ^ - "^{n] - Uj) - n\ (^^^ + ■ (4-28) 

The choice iXj = —{fij — l)/2 is thus the only one such that the argument of 
the exponential in (14.271) is a real number. It is also the only choice such that 
'Pk{x'i M, n,q + inpi) and ^niy, M, n,q + infl + i/tn — in) are complex conjugates of each 
other. Introducing the Bethe eigenfunctions [11] 

x;M,n,q--n + -\ , (4.29) 

one finds for fij = —{rij — l)/2 



1 f 



M=l ^ ' n€D„^M 

M 

ij{x- M, n, q) ij{y;M,n, q) JJ e-*"^'?^'^^ K'""^) . (4.30) 
This is the result (13. 8 p of Theorem |2l 
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5. The special case x = y = Q 

In the special case a; = ?/ = 0, Eq. fl4.30p simplifies. One has 

' {qj + iK(/Xj + r(a))) - {q^ + ift:(/Xfc + r{h))) + ik 
{qj + + r(a))) - (g^ + in^fik + r{b))) - in 
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M 



E n n 

aeS^{n)j<k aesij 
A/ 



(5.1) 



E n n ((gj -gfc) + ift(/ii -/ifc + r(a) -r(6) + sign(a H«) H^)))) 



cTeS'{n)j<k aen^ 

b 6 Qfc 



M rij uk y 

n n n ( - + i'^i/^j - f^k+r-s-r 

i<kr=ls=l ^ 



As in the previous section, one can sum instead over all permutations a G S'„ in case the 
product in the numerator becomes a product over all a < b: the permutations which do 
not belong to S'^{n) give zero contribution. We use the notation C,a = Ij + i'^ifJ'j + r(a)) 
for a & Qj. Since the signature of the permutation a can be written 

fj ^a-^b + iti sign((T"^(a) - cr^^{b)) 



sign (^0- J 



n 

a<b 



^a(a) - ^a(fe) + sign(a - 6) 



(5.2) 



Eq. ( 15. ip rewrites as 

M \ E sign(cr) n Uaia) - ^a(fe) + i/t sign(a - 6) 



nj(nj — l) 



-IK) 



nn 



L J. J. I M «j rife X X 
lr-=l / ]^ J] J] ^(g^. + 

j<k r=l s=l ^ ^ 



(5.3) 



One has {e.g. see [19], Lemma 1 for a proof) 

Lemma 3. Let f{a,b) be arbitrary complex coefficients. Then 

n n 

sign((T) Yl (^a(a) - ^a(b) + f{a, b)) = n\ Jli^a - 6) • 

adSn a<b a<b 

Using this lemma, (15. 3p rewrites 



(5.4) 



, / T-r -i-r 1 



n\(—iKi 



nn 



n (u - &) 

a<b 



n;^ nnn 



■I / M nfc 

n 

j<k 1 

{(Ij - (Ik) + i/t(/ii - lik + r - s) 



ir-=i •/ Y[]\]\[{qJ-qk) + H^^J-^^k + r-s-l)) 

i<k r=l s=l ^ ' 



\j=l ■' / j<ck r= 



,1 fli V - Qk) + i/t(/ij - /ife + r - s - 1) 
Then using (see [H], Eq. (B.55-B.58) and Eq. (34)) 



(5.5) 



M rij rik 

nnn 

j<k r=l s=l 



{qj - Qk) + i/t(- 



+ ^ + r - s 



(gj - gfe) + i«:(-T + f + r - s - 1) 
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one finds 



M 

n 

det 



- (Ik) - f K- - nk) 



mrij 



{<lj - Qk) + f{nj + nk) 



j,k=l,...,M 



M-- 



det 



M\{2ti 



I A/ 



dqi ... dq^ 



Q-tnjq^ + !%-(n^^-n,) 

-KQj - Qk) + ^{nj +nk) 



M 



j,k=l,...,M 



(5.6) 



(5.7) 



In the replica computations for tlie KPZ equation witli sliarp wedge initial data, Eq. 
f l5.7p can be used to arrive directly at the generating function of the height fluctuations, 
see Eq. (2.17)-(2.19) of [I9]. 



6. Conclusions 



Exact formulas for the transition probability of the one- dimensional asymmetric simple 
exclusion process, a non-equilibrium exactly solvable model, have been derived a few 
years ago [321 [331 [31] . Subsequently the method was adapted to obtain an exact formula 
for the propagator of the quantum (5-Bose gas with repulsive interaction [26]. Here we 
analytically continued this formula to the case of attractive interaction. 

An advantage of our approach, compared to the usual Bethe ansatz, is that the 
question of the completeness of the Bethe eigenfunctions can be completely avoided. 
In fact, such kind of exact expressions for the propagator can be used to prove the 
completeness of the Bethe ansatz, at least on the infinite line. In principle, it might 
also be possible to use the same kind of approach for the case of periodic boundary 
conditions, see e.g. [331 13^- 

The exact expression for the propagator of the repulsive (5-Bose gas, derived in 
[26] and used here, bears some formal similarity with the coordinate Bethe ansatz, 
which is the original ansatz introduced by Bethe to diagonalize the Hamiltonian of 
Heisenberg spin chain. Since then, other descriptions of eigenstates have been developed, 
in particular algebraic Bethe ansatz, which makes clearer the structures underlying the 
quantum integrability of such type of models. It would be of interest to understand 
whether it is possible to write down the propagator using an approach closer to the 
algebraic Bethe ansatz. 
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Appendix A. Analyticity in the coupling of the propagator of the 5-Bose gas 

In the case of a standard Schrodinger operator of the form — A + W one can use Kato's 
theory to establish that e"*^"^"'"''*^) is analytic in A. The 5-potential corresponds to 
a boundary condition and we are not aware of a functional analytic argument for the 
holomorphic dependence on k. Instead, we will use the Feynman-Kac representation. 

Proposition 1. For fixed x,y G M", t > 0, the function k, i— )■ (xle^^^^ly) is holomorphic 
on C. 



Proof. By the Feynman-Kac formula one has the representation 

(x|e-*^1y) = E,,y (e^'^^W) p,{x - y) . (A.l) 

Here pt{x — y) is the Brownian motion transition kernel. The expectation K^^y is over 
the standard Brownian bridge, b{t), starting at x and ending at y at time t. Let Lj^k{t), 
j < k, be the local time at for — bj{s),0 < s < t}, i.e. 

L^,k{t)= I ds5{h^{s)-h{s)) . (A.2) 

Then 

n 

X{t) = Y,L,,u{t). (A.3) 

j<k 

Since X[t) > 0, we bound as 

n 

E.,. (e^''^^*^ < E.,. (e^l'^l^^*)) < E,,, (el'^l"^^^.^ W) . (A.4) 

j<k 

The second inequality follows from the convexity of the exponential function. Lj^k{t) is 
the local time at for a one-dimensional Brownian bridge starting at xj — Xk and ending 
at yj — yk at time t. The distribution of Lj^k{t) has a Gaussian decay at infinity. Hence 
for any \k\ the right hand side in ( 1A.4I) is bounded. This proves analyticity on C. □ 
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